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We study the geometrical corrections to the simple Proximity Force Approximation for the non- 
retarded Casimir force. We present analytical results for the force between objects of various shapes 
and substrates, and between pairs of objects. We compare the results to those from more exact nu- 
merical calculations. We treat spheres, spheroids, cylinders, cubes, cones, and wings; the analytical 
PFA results together with the geometrical correction factors are summarized in a table. 
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I. INTRODUCTION 

J. D. van der Waals found empirically in 1873 that there is an attractive force between non-polar atoms. It took 
^""i a long time, until 1930, before there was an explanation to this force. London [11] gave the explanation in terms of 

^_^ fluctuations in the electron density within the atoms (fluctuating dipoles). In an alternative description ^2] one may, 
r^ , instead of discussing the particles, focus on the electromagnetic fields. The force may be expressed as a result of 
rS i' changes in the zero-point energy of the electromagnetic normal modes of the system. There are modes associated 
•j , with the atoms and modes associated with the vacuum. Casimir studied a more pure and idealized system consisting 
^ • of two perfectly reflecting metal plates. In this geometry there are no modes associated with the plates themselves; 
2 ' there are only vacuum modes. The presence of the plates changes the vacuum modes and their zero-point energy. He 
^ |: published his findings in a classical paper '31 in 1948, the Casimir force was born, 60 years ago. Casimir's paper is one 
of the most important papers in the history of physics since it demonstrates that the boundary conditions of a system 
may change its zero-point energy and hence its properties. In the case of objects made from real materials both the 
modes associated with the objects and the vacuum modes contribute to the force 2] . The first type dominates in the 
van der Waals region which is for smaller separations. For larger separations between the objects the vacuum modes 
dominate and the result is the Casimir force. In the idealized case of perfect metals there is no van der Waals range. 



o^: 



> 



\f^ ' In 1997 Lamoreaux [4] performed the first modern high-precision measurements of the Casimir force. The accuracy 
was good enough to make the direct comparison with theory feasible. This spurred a burst of renewed interest in the 
Casimir effect. 

The dispersion forces (van der Waals and Casimir forces) decrease in size with separation much faster than the 
QQ ' Coulomb force and the gravitational force do. This means that for macroscopic objects at macroscopic separations 
f~^ , the dispersion forces are very weak compared to other forces. However, in the micrometer and nanometer ranges 
they often dominate. Also in biological systems they are very important; the Coulomb forces are screened by the 
abundance of mobile ions; the dispersion forces are not screened. 

In the rapidly emerging field of nano-technology the dispersion forces have come to play an important role and one 
5— i ' tries to exploit these forces in nano-mechanical devices. In the design of such devices it is important to develop precise 
calculations of the Casimir forces between objects with different geometries and in different configurations. It may be 
complicated to calculate the force between objects of general shapes. The van der Waals and Casimir forces between 
half spaces is, however, not very difficult to calculate. For finite objects at large distances one may use multipolar 
expansions [5|. However, for small separations one needs to keep more and more terms in the expansion the smaller 
the distance; one reaches a limit when the method is no longer feasible to use. 

In this work we study an approximation that comes handy in this situation. The Proximity Force Approximation 
(PFA). It was first used already in 1934 [6] in connection with coagulation of aerosols. It is a very powerful and widely 
used approximation for the interaction at short distances between two objects. Lamoreaux' experiment, discussed 
above, is interpreted within this approximation. It is difficult to make a strict estimate of how good the approximation 
is but it has gained a wide-spread acceptance in recent years. 

The basic idea of the approximation is that the interaction potential between the objects is an average interaction 
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energy between parallel planar interfaces, 

V{z)^ I dSEpiw), (1) 



where Ep (w) is the interaction energy per unit area between planar interfaces a distance w apart. The variable z 
is the closest distance between the two objects. The surface S is not uniquely defined and the choice of S is in the 
general case not a trivial choice to make. In the examples we discuss here the choice is more obvious. For objects 
above a substrate we may choose S to be that part of a flat surface, parallel to the substrate, that is covered by the 
projection of the object; w is the distance at dS between the substrate and the object along the normal to S. In the 
case of two interacting objects we suggest that S is that part of the planar surface, perpendicular to the resulting 
force, where the projections of the two objects overlap; w is the distance at dS between the two objects along the 
normal to S. 

Deviations from the PFA results have been observed in experiments dealing with spheres. Q There have been 
suggestions that the deviations are due to the thin coatings of the objects. However due to limitations in the 
experimental and theoretical parameters one can not say conclusively that the PFA calculations for coated plates give 
better predictions than PFA for solid and infinitely thick plates. 

Some attempts have been made to show its validity in different cases. However, the lack of a general and rigorous 
proof still remains, as well as estimates of the level of accuracy at all distances. In the present work we follow an 
alternative procedure where we perform an exhaustive study of the PFA results for a wide range of geometries and 
configurations, and compare them with the most exact results available in the literature. 

A general method for the calculation of van der Waals forces is shown in Ref. [g]. It is based on the solutions 
of a surface integral equation for the electric potential. It is fast and robust, and allows for calculations with great 
accuracy. We use either this method or a multi-polar expansion method to compare with the PFA results. The multi- 
polar expansion method we use in the case of spheres and spheroids where the objects have simple symmetries. There 
are two versions of the integral equation method, one three-dimensional (3D) and one two-dimensional (2D). The 3D 
version is the more general of the two. In some situations one may use the 2D version which is faster converging. 
An example where this is possible is a finite cylinder above and parallel to a substrate. The force is in this case to a 
good approximation equal to the length of the cylinder times the force per unit length for an infinite cylinder above 
a substrate. In the infinite cylinder case we have cylinder symmetry and a 2D problem. The 3D and 2D versions give 
the same result for small separations but give different results for large. When they deviate it is the results from the 
3D version that is correct. 

We will in the sections that follow go through different examples and compare the results to those from more 
elaborate calculations. We use extensions to the traditional PFA as discussed in Sec. |TT1 In Sec. Illll we show how the 
calculation of the interaction energy per unit area in planar structures are performed. These results are used in our 
calculations for more general structures. In Sec. IIVI we give a general expression for the interaction energy for half 
spaces, cyhnders and spheres. Sec. |V]is devoted to spherical objects, Sec. |Vl]to spheroids. Sec. I VIII to cylinders. Sec. 
IVIIII to cubes. Sec. llXl to cones, and Sec. |X]to wings. Finally, summary and conclusions are found in Sec. IXII The 
main analytical results are summarized in Table U of Sec. IXII 



II. THE PROXIMITY FORCE APPROXIMATION IN THIS WORK 

Traditionally in the PFA of the interaction between two objects one only takes into consideration the surface of 
each object that is facing the other object. The function Ep (w), in Eq. ([T]), is the energy per unit area for two half 
spaces, made up from the materials of the two objects, separated by the distance w, see Fig. [TJa). If the objects 
are immersed in an ambient medium the gap between the two half spaces is also filled with this medium. With this 
treatment the backsides of the objects have no effects at all. One gets, e.g., the same result for two spheres as for 
two halfspheres. We know that the normal modes contributing to the interaction may extend through the objects 
and continue on the other side. Thus, the backside may have important effects on the results. One may extend 
the treatment by using Ep (w) from a planar structure with four interfaces instead of two, according to Fig. [IKb). 
The distance between the interfaces surrounding the object material is the local thickness of the object. We will not 
consider this in the present work, but address another limitation that is more severe. Two coated spheres, e.g., are 
in PFA treated as solid spheres. We will extend the treatment by considering a planar structure with four interfaces, 
like in Fig. [Ijc). In the most general case the two coatings are of different materials and of different thickness. We 
let Ep {w, 5) be the energy per unit area for two coated half spaces, made up from the materials of the two objects, 
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FIG. 2: Coated oblate spheroid above a substrate used to 



FIG. 1: (a) Planar structure entering Ep (w) in the tradi- 
tional PFA. (b) Planar structure used in an extension of illustrate the parameters discussed in the text. 
PFA to take the back sides of the objects into account, (c) 
Planar structure used in the calculation of Ep {w, S) in our 
treatment of coated objects. 



separated the distance w. The thickness of the coating is denoted by 6. We may express the potential as 



2 + A 



V{z)^ / dSEp {w, 6) 



dw —— Ep (w, 6) , 
aw 

g{w-z) 



and find the force as 



(2) 



F{z) = -^ = -g{^)Ep{z + ^,5)+g{Q)Ep{z,5)- 



dg{w~ z) 

dw ; Ep [W, 0) . 

dz 



(3) 



In many cases g (A) vanishes, like in the illustrating example in Fig. [S] Here we have a coated spheroid above a 
substrate. The surface S is defined through the projection of the object onto the substrate. When g (A) vanishes we 
have 



F{z)=giO)Ep{zJ) 



1 



z+A 



g{0)Ep{z,6) 



fPS 

dw--—Ep (w, S) 
dw^ 



(4) 



where the first part is what one usually means with PFA. The remaining part, within brackets, is a correction factor 
depending on the geometry. This factor is often dropped without any motivation at all or with the argument that 
the resulting error is of the same order of magnitude as the error in PFA itself. Often Ep (w, S) has a known power 
dependence in w over the whole integration interval. Then we may often find an analytical expression for the correction 
factor. In Fig. [5] we study a coated oblate spheroid above an uncoated substrate. The lowest horizontal line is the 
substrate boundary. The thicker part of this line indicates the surface S, which is defined by the projection of the 
spheroid on the substrate. The integration variable w varies from z, the closest distance, to z + A. The region 
between w and w -\- dw picks out a surface on the object whose projection on the substrate is a ring of area dS. This 
surface contributes, in the PFA spirit, to the energy with dS times the energy per unit area, Ep {w, 6), of the planar 
configuration defined by the substrate and the film of thickness 5, indicated in the figure. If the spheroid is not an 
empty shell, the half space above the film should be filled by the material of the spheroid. 

To summarize, in this work we extend the traditional PFA in two ways; we retain the correction factor of Eq. ([¥]); 
we extend PFA to include the finite coat thickness of coated objects. The effects of these extensions are demonstrated 
and comparisons are made to numerical results from more accurate calculations. Throughout the text we refer to 
results from traditional PFA as PFA results and to extended or corrected PFA results as full PFA results. 

Before we proceed with the various objects we derive, in next section, the interaction energy in those planar 
structures we need for the PFA calculations. 



III. PLANAR STRUCTURES 

A general expression for the interaction energy per unit area in a planar system where only one of the distances, 
w, between neighboring interfaces is allowed to vary is [2] 

EpH^'j7 tlS^ {In [/ (k, *^, w)] - In [/ (k, lu;, oo)]} 

-^ . .,.. ,. (5) 



h r duj_ r d k 1 

2 J 27r J (2^ 



/(k,ia;,oo) 



where the reference energy is set to when w is infinite. The variable k is the two dimensional wave vector in the plane 
of the interfaces. The function / (k, WjUi) is the function in the condition for having electromagnetic normal modes 
in the planar system, 

/(k,c^,«;)=0. (6) 

In the general case there are two of these functions, one for TE modes and one for TM modes. Here we limit the 
treatment to the non-retarded limit. One of the effects from neglecting retardation is that the TE modes are absent. 
In the case of a coated object above a coated substrate or of two coated objects we need four interfaces between five 
regions. For a structure of the type 1|2|3|4|5 where medium 3 has the variable thickness w we have 

f (k, CJ, oo) ' 1 + e-2fed2 j,^2r23 + e-2'='i4,^3^r45 + e-2'='i2e-2fcd4,^^2r23r34r45 ' 

where /c = |k| and 

£j- ju;) - £, {iu) 
"^'^-e.H + e.H- ^^^ 

The function Si (uj) and di are the dielectric function and thickness, respectively, of region i. In the case of a coated 
object above a substrate or of two objects where only one is coated we need three interfaces between four regions. 
For a structure of the type 1|3|4|5 where medium 3 has the variable thickness w we let ri2 = and get 

/(k,^,w) ^ ^ ^ ^_2fc«, n3?-34 + e-^'"^->ri3r45 ,^. 

f{k,uj,oo) 1 + 6-2^^47.34^45 ■ ^ ' 

In the case of an uncoated object above an uncoated substrate or of two uncoated objects we need two interfaces 
between three regions. For a structure of the type 1|3|5 where medium 3 has the variable thickness w we let r^^ — 
and get 

f(\i.,LU,w) oj,,, 

i;, \ = 1 + e-2^-ri3r35. 10 

/(k,w,cx)) 

For this last structure the integration over momentum in Eq. ([5]) may be performed and results in an infinite series. 
Ep {w) = ^ J dcuj dkkln [1 - e-2fe»^^3^53] ^ _^ J du; J dkkln [l - e-'^n^rsa] 



"327r^u,i J "''^ 2^ P - ^327r'^w-^ 2^ ' 

-oo 1=1 1=1 



(11) 



where the characteristic dielectric integrals (uji) are 

oo 

(uJi) ^ duj 



£l (uj) - £3 (w) £5 (w) - £3 (w)"' 



£1 {uj) + £3 (uj) £5 (u;) + £3 (uj) 



(12) 



This result is obtained after a variable substitution, series expansion of the logarithm, followed by the integration 
over momentum. 
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FIG. 3: The correction factor CwriC,, x) for a set of C, values. 

IV. GENERAL EXPRESSION FOR HALF SPACES, CYLINDERS AND SPHERES 



For solid objects in neglect of the geometrical correction in Eq. ([4]) one may find the following general expression 
for half spaces, cylinders and spheres: [2, 0] 



E{z) 



327r2z(i+"/2) 



T {1 + 71/2) 



2ttRiR2 
R1+R2 



1 — n/2 00 / , 

1=1 



(13) 



where n = for spheres, n — 1 for cylinders, and n = 2 for half spaces. The result is for the interaction energy in 
case of spheres, the interaction energy per unit length for cylinders and the interaction energy per unit area for half 
spaces. The variable z is the closest distance between the objects. Thus, for two half spaces we have 



for two cylinders of radii i?i and i?2 



E{z) 



r(3/2) 



327r2z2 ^ P 



±^4=M^): 



2ttRiR2 

Rl + i?2 



1/2 00 



1=1 



P 



i:¥=r(3/2),/lssiE,(,,, 



and for two spheres of radii Ri and R2 the result is 



E{z) 



327T^Z 



2ttRiR2 



Rl + i?2 



E 

1=1 



^3 



2ttRiR2Z 
Rl + R2 



Rl + i?2 



Ej, (z) . 



(14) 



(15) 



(16) 



To get the result for a cylinder of radius R above a substrate we let R2 go to infinity and replace Ri with R in Eq. 
(fTSll. This results in 



E{z)^- 



^ F (3/2) V2^ E "ir = ^ (3/2) V2^zEp {z) . 



327r2z3/2 



1=1 



P 



(17) 



To get the result for a sphere of radius R above a substrate we let R2 go to infinity and replace Ri with R in Eq. 
[TBI). This results in 



E{z) 



^ 27rRf2^ = 27rRzEp{z). 
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1=1 
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(18) 



This is as far as we get with the general formula. Now, we continue in next and the following sections with the 
geometrical corrections. We start with spherical objects. 
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FIG. 4: R times the non-retarded force on a gold sphere 
above a gold substrate as function of z/ R for the multipolar 
result, the PFA result and the geometry corrected or full 
PFA result. 
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FIG. 5: The non-retarded force on a gold sphere above a 
gold substrate relative the PFA result as function of z/R for 
the multipolar result and the full PFA result. 



SPHERICAL OBJECTS 



Spherical objects are often used in experiments. The advantage is that one avoids the problem of alignment. This is 
the case both for a sphere above a substrate and for two interacting spheres. In Lamoreaux's [^ classical measurement 
of the force between two gold plates one actually measured the force between a sphere and a planar surface. 



A. Sphere-substrate interaction 

For a sphere above a substrate the parameters entering Eq. ^ are 

A = i?; g{x)^2T:{R-x)- g (0) = 27ri?; ^^ = -2^, 



and this results in 



F (z) = 2ttREp {z, 5) 



z+R 



1- 



REp (z, ,5) 



dwEp (w, S) 



(19) 



(20) 



where we have included the possibility for the sphere and/or the substrate to have a coating of thickness 6. If Ep {w, S) 
varies as 



Ep {w, S) = -C/w^ 
in the whole integration interval we find the correction factor, Corr{C,, x) on analytical form 



Corr (C, x) 



1- 



x^/{l + x)^-']/{C-l): Cy^l 






l-x\n[{l + x)/x\; C = l 



(21) 



(22) 



where x = z/R. The correction factor is illustrated for some C, values in Fig. [3l Eq. ([20| is valid both in the 
non-retarded and retarded separation regions. The PFA is only good for small x-values but if R is big enough one 
may still be in the retarded region. For uncoated sphere and substrate ^ = 2 (C = 3) in the non-retarded (retarded) 
region. If either of the sphere and substrate or both are coated the separation dependence is in general more complex. 
If a metallic coating is thin enough there is a separation range where Q = 5/2. 

With the proper scaling one may produce universal figures, i.e., figures that are identical for all values of R. From 
Eqs. © and follows that Ep{w,6) = f {w/R,S/R)/w^. From Eq. ^ follows that if one plots i? x F as a 
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FIG. 6: R times the non-retarded force on a spherical gold shell above a gold substrate as function of z/R. The curves are for 
the multipolar result (thin solid curve with open circles), the PFA result for a solid sphere (thick solid curve), the full PFA 
result for a solid sphere (dashed curve), the two-dimensional PFA result (dotted curve), the PFA result for the shell (solid 
circles) and the full PFA for the spherical shell (triangles). 



function of z/B. the resuhs look identical as long as the coat thickness is the same fraction of B. This means that 
there is no need to repeat the calculations for different sphere radii. In Fig. [4]we show B times the non-retarded force 
on a gold sphere above a gold substrate as function of normalized separation z/B. Throughout we use the dielectric 
function of gold given in Ref. [lOl]. The differences between the curves are not clearly seen in a figure like this. In 
Fig. [5] we show the result from the multipolar calculation and the geometry corrected or full PFA result relative the 
PFA result. We see that the full PFA means a substantial improvement from PFA. In Fig. [6] we show B times the 
non-retarded force on a spherical gold shell above a gold substrate as function of normalized separation z/B. The 
thick solid curve and dashed curve are the PFA and full PFA curves, respectively, for a solid gold sphere and are 
the same as in Fig. 2] the solid curve with open circles is the result from the multipolar calculation for a shell of a 
thickness of one percent of the radius; the filled circles are our extended PFA result for the gold shell; the solid curve 
with triangles is our extended and full PFA result for the gold shell. We note that the full result, from the multipolar 
calculation, follows the PFA result for a solid gold sphere for distances smaller than approximately the coat thickness. 
Then for larger separations there is a region where it follows rather closely the 2D PFA result obtained by replacing 
the gold film with a 2D (two dimensional) metallic sheet with 2D electron density given as the projection of the 3D 
electron density of the film. These results are derived in Refs. |lll.[l2|. 



Bx F ::^ 2t:B^Ep {z) w O.lbbQ^nfi^e'^/mey/sjB {z/Rf^, 
where n is the conduction electron density of gold. 



(23) 



B. Sphere-sphere interaction 

For two spheres of equal size the parameters entering Eq. ([4]) are 



A = 2i?; g{x)=n{B-x/2)- g(Q) = ^B- 






^V2, 



and this results in 



F (z) = irBEp {z, 6) 



Z+2R 



dwEp (ui, 5) 



(24) 



(25) 



2BEp {z, 6) 
If Ep (w, 6) varies as in Eq. (|2ip in the whole integration interval we again find the correction factor, Corr{C,, x), on 
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FIG. 7: 8R times the non-retarded force between two gold 
spheres as function of z/2R for the multipolar result, the 
PFA result and the full PFA result. 



FIG. 8: The non-retarded force between two equal size gold 
spheres relative the PFA result as function of z/2R for the 
multipolar result and the full PFA result. 



analytical form and it is identical to the expression in Eq. (P^ but now x = z/2R. In Fig. [7] we show 8R times the 
non-retarded force between two gold spheres of equal size as function of normalized separation z/2R. Plotted in this 
way the PFA and full PFA curves are identical to the corresponding ones in Fig. 31 the result from the multipolar 
calculation is not. In Fig. [8] we show the result from the multipolar calculation and the full PFA result relative the 
PFA result. We find that the agreement is not as good as in the case of a sphere above a substrate. 



VI. SPHEROIDS 

A spheroid has an axis of symmetry. We let this axis be perpendicular to the substrate in case of a spheroid above 
a substrate; in the case of two spheroids we let both symmetry axes coincide with the line joining the centers of the 
spheroids. We let y be the variable along the symmetry axis and r the variable perpendicular to the symmetry axis. 
Let A (B) be the larger (smaller) of the two semiaxes and let 



7 = A/B > 1. 
There are two types of spheroid, oblate defined by the equation 

2 9 

r y 

and prolate defined by 



(26) 



(27) 



r 



A2 



= 1. 



(28) 



The oblate spheroid varies between a sphere of radius A when B is equal to A and a disk of radius A when B is equal 
to 0; the prolate spheroid varies between a sphere of radius A when B is equal to A and a pin of length A when B is 
equal to 0. 

A. Spheroid-substrate interaction 

We start with the oblate spheroids. For an oblate spheroid above a substrate the parameters entering Eq. ([4]) are 



A = B; g{x)^-2^-i\x-B)- g (0) ^ 27rj^B; ^ = -2^^\ 
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(29) 
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FIG. 9: The force on a gold spheroid above a gold substrate. 
The thick solid curve and dashed curve are for the PFA and 
full PFA, respectively. The upper (lower) solid curve is for 
a prolate (oblate) spheroid with 7 = 4.0 (7 = 1.4). The 
left (right) and lower (upper) axes are for oblate (prolate) 
spheroids. 
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FIG. 10: The force between two gold spheroids. The thick 
solid curve and dashed curve are for the PFA and full PFA, 
respectively. The upper (lower) solid curve is for a prolate 
(oblate) spheroid with 7 = 4.0 (7 = 1.4). The left (right) 
and lower (upper) axes are for oblate (prolate) spheroids. 



and this results in 



Thus, 



F (z) = 2TrB"/^Ep (z, S) 



1- 
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z+B 



BEp {z, 5) 



dwEp {w, S) 



F (z) = 27rB-f^Ep {z, S) Corr (C, x) ; x = z/B, 

if Ep (z, 6) has a power law dependence. The correction factor, Corr {(^, x), is the same as in Eq. (|22 
To get the results for prolate spheroids we just replace B with A and 7 with 7^^. The result is 



F (z) = 2ttAj-^Ep (z, d) 



z+A 



AEp (z, 6) 



and if Ep (z, S) has a power law dependence it is 

F (z) = 2TTA-f-'^Ep (2, S) Corr (C, x) 



dwEp {w,S) 



z/A. 



(30) 



(31) 



(32) 



(33) 



The correction factor, Corr {C,x), is the same as in Eq. ((22|) . 

If one for an oblate (prolate) spheroid plots B-f^^ x F (^7^ x F) as a function of z/B {z/A) and lets the thickness 
of the coating be in units of B {A) one obtains the same universal PFA and full PFA curves as we previously found for 
a sphere above a substrate or for two spheres. Fig. [5] shows the results for an oblate (prolate) spheroid with 7 = 1.4 

(7 = 4.0) above a gold substrate. An alternative way to plot the curves is to plot R x F as function of z / R, where 

R is the radius of curvature at the point of closest contact. Then the PFA result is universal but the full PFA result 
is different for the oblate and prolate cases and depends on the 7 value. The radius of curvature is Bj^ (^7"^) for 
oblate (prolate) spheroids. 



B. Spheroid-spheroid interaction 



In the case two equal spheroids aligned along the common symmetry axis we find in the oblate case 

Z+2B 



F (z) = TTBj^Ep (z, S) 



1 



1 



2BEp (z, S) 



dwEp {w, 5) 



(34) 



and if Ep (z, S) has a power law dependence it is 

F (z) = TrBj'^Ep {z, 5) Corr (C, x) ; x^ z/2B. 

The correction factor, Corr (Cx), is the same as in Eq. 
For two prolate spheroids we obtain 
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(35) 



F (z) = TTA-f-^Ep {z, S) 



-2A 



1- 



2AEp (z, 5) 



dwEp (w, 5) 



and 



F (z) = nA'y-'^Ep (z, ,5) Corr (C, x) ; x^ z/2A, 



(36) 



(37) 



if Ep (z, S) has a power law dependence. The correction factor, Corr {(, x), is the same as in Eq. ([22|) . If one for two 
oblate (prolate) spheroids plots 8Bj~'^ x F (8^7^ x F) as a function of z/2B {z/2A) and lets the thickness of the 
coating be in units of 2B {2A) one obtains the same universal PFA and full PFA curves as we previously found for a 
sphere above a substrate, for two spheres and for a spheroid above a substrate. In Fig. [10] we show the force between 
two gold spheroids, two prolate with 7 = 4.0 and two oblate with 7 = 1.4. An alternative way to plot the curves is to 
plot 8R X F as function of z/2R, where R is the radius of curvature at the point of closest contact. Then the PFA 
result is universal but the full PFA result is different for the oblate and prolate cases and depends on the 7 value. 

VII. CYLINDERS 

The problem of alignment in measurements of the dispersion forces is partly avoided by using cylinders. It is fully 
avoided in case of crossed cylinders jlSj 



A. Cylinder-substrate interaction 

Here, we consider a cylinder of radius R and length L above a substrate. The closest distance is z, as before. This 
is a geometry where once again g (A) of Eq. ([3]) vanishes. However, here one runs into another problem, viz., that 
g (0) diverges. We need to make an alternative derivation. 



y (z) = / dSEp (w) = 2LJ dyEp [z + R- ^. 



2_y2 



2L / dw 





+rt.-w 



w = z + R- yjR^ - y2 



y/R:^-{z+R-w) 



--Ep (w) = \w ^ u + z\ — 2L J du 



(38) 



R-u 



and 



F{z) 



dV{z) 
dz 



R 

'2L / du- 



R-u 



^ ^/R^-(R^u) 



dEp {u + z) 



=Ep (u + z) , 



dz 



Vi?2-(i?-u) 

We limit the treatment to power law dependence in the non-retarded limit and obtain 



F(z)= LEp{z){i^/A)./2^^z 

[•iyfTT / i)^ 2tI / X[x+l)(x + 2.f ^ X'^ +2X 



(39) 



(40) 



F{z)R^/L= Fp (x) (30F/4) ^2^ 

Q(x+lf tan~'^{ll ^x'^+2x)+Z-jr[x+l)'^+2^x'^+2x{2x'^+Ax+3,) 
(3V^/4)-/27i7a;(a;+l)(a:+2)^VaJ^+2a; 



(41) 
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FIG. 11: The force on a gold cylinder of radius R and length L above a gold substrate. The thick solid curve and dashed curve 
are for the PFA and full PFA, respectively. The thin solid curve is the full result from calculations based on integral equation 
techniques. 

where we have written the results on universal form. The first factor is the PFA result and the second the correction 
factor from the full geometry. The variable x = z/R. The result is valid under the assumption that E,p (z) ex z^^. Fig. 
[TT] shows the result (thin solid curve) for a gold cylinder above a gold substrate as calculated with integral equation 
techniques [8|. We have performed calculations for an infinite cylinder and for finite ones with L = R and L = 2R. 
All three curves fall within the curve thickness on a plot like this. The result fronr the infinite cylinder calculation 
we limited to z < R otherwise we would have gotten deviations. The thick solid (thin dashed) curve is the PFA (full 
PFA) resuh. 



B. Cylinder-cylinder interaction 



The result for the force between two identical parallel cylinders is easily obtained from the results in the previous 
subsection. It is 



F (z) 8R'^/L = Ep (x) (30F/4) ^2^2? 

6(2; + r)^ tan"i(l/V2;2+2a;)+3ir(£c+l)^+2Va;^+2a;(2£i;^+4x+3) 

where now the variable x = z/2R. Note that the right hand side of this equation is identical to the one in Eq. (|4T 



(42) 



C. Standing cylinder 



For a cylinder standing upright above a substrate, the proximity force approximation gives the potential as the 
bottom area of the cylinder times the interaction energy per unit area between two halfspaces. The full PFA gives 
the same result. For a cylinder of length L with circular cross section of radius R, we have in the non- retarded case 
V (z) ^ TiR'^Ep (z) and 



The universal expression is 



F{z) = -^^ = 27rR^-Epiz). 
dz z 



RF = 27rEp {x)/x, 



(43) 



(44) 



where x = z/R. 
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FIG. 12: The force on a gold cube of side length L above a gold substrate. The dotted curves and dashed curves are for the 
PFA and full PFA, respectively. The thin solid curves is the full result from calculations based on integral equation techniques. 

VIII. CUBE-SUBSTRATE INTERACTION 

A. Straight cube 

For a cube standing upright above a substrate, the proximity force approximation gives the potential as the bottom 
area of the cube times the interaction energy per unit area between two halfspaces. The full PFA gives the same 
result. For a cube of side length L we have in the non-retarded case V (z) = L^Ep (z) and 



The universal expression is 



where x ~ z/L. 



dViz) ,1 

Fiz) = -^^2L--Epiz) 



LF = 2Ep (x)/x, 



(45) 



(46) 



B. Tilted cube 



Here we tilt the cube the angle a while keeping one of the edges parallel to the substrate. The potential is 



and the force 



V {z) ~ j dSEp(w) — J dxLEp (z + xtana) + J dxLEp(z + x cot a) 
so 

z-^L sin a z+L cos a 

= Lcota J dwEplw) + L tan a J dwEp{w), 

z z 



F{z)^- rf^ ^ _^ [(,Q^ ^^^ (z + L sin a) - cot aEp (z) 
+ tan aEp {z + L cos a) — tan aEp (z)] 
= sina^osa ['^o^^ ^^P (^ + L sma) + siu^ aEp {z + L cos a) - Ep (z) 



The universal expression in the non-retarded limit is 

1 



LF 



sm a cos a 



-Ep [x) 



1 - cos a 1 



sma 



sin a I 1 



/ cos a \ 



(47) 



(48) 



(49) 



13 



> 






10' 
10' 
10" 

io-> 

10"^ 
10-^ 
10-' 













V 


— 


i / - 














\ 


L 


a / 
















\ 




*^^5^^^ PFA 














\^ 


^ ^^^^^^ \ 














^t 


^^ ^^^^^^^ 


\ 














^ 




^ 










- 


\ 








^ 






- 


^ 


^ 












.^ 






\ 


■^ 


\ 


\ 




Full PFA - 



































m 



10" 



10^ 



z/L 



FIG. 13: The force on a gold cone of height L above a gold substrate. Plotted in this way the PFA and full PFA results (thick 
solid and dashed curve, respectively) produce universal curves independent of scale and aperture angle (2a). 

where x = z/L. The first factor is the PFA result and the second the correction factor due to the full geometry. 
Note that the x dependence of the PFA result changes abruptly when the cube is tilted. In Fig. [T^] we show the 
results for a gold cube above a gold substrate. The results are for a straight cube and for a cube tilted by 7r/9 and 
7r/4, respectively. The dotted curves are the PFA result and the dashed curves the full PFA result. The results from 
the integral equation method are represented by thin solid curves. For the straight cube there are no corrections to 
PFA. For large separations the force for both straight and tilted cubes varies as F ^ L^/z^; for small separations the 
force on the straight cube varies as F ~ I? j z^\ for small separations the force on a tilted cube varies as F ~ Ljz^. 
This means that the force at large separations is proportional to the volume of the object; at small separations it is 
proportional to the projected area for the straight cube and to the length of the closest edge for the tilted cube. 



IX. CONE-SUBSTRATE INTERACTION 
A. Cone with a pointed tip 



Let the cone have circular bottom of radius R and be of height L. It has its point towards the substrate at the 
distance z. In this case Eqs. ^ and ^ can not be used. Instead, we write 



Vih)^ f dSEp (w) = I dr2TrrEp {h + rL/R) , 



(50) 



and if Eq. pT|) holds we have 



F{z) 



dV{z) _ 



dz 



= -Jdr 



27rrCn 



{z+rL/Ry 



2TTCnR 



?^Idn 



^ ' ' 



-27rCni?2 [dr^—!- 

-> iz-\-rL 



{z+rL) 



ITT 



2'kCR' 



1 

n-1 



{z/L + n)/{z/L+lT-l {z/LY 



(51) 



n^l 



[ 1 + In [{z/L) /{z/L + 1)] - {z/L) /{z/L + 1); n = 1. 
In our case we have n = 2 and may write on universal form 

F{z)L = 2TrEp {x) x^ tan^ (a) / x (1 + x)^ ; x = z/L. 



(52) 
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FIG. 14: The force on a gold cone of height L with a spher- 
ical tip with radius of curvature R above a gold substrate. 
The results are for L/R = 8. The thin solid (dotted) curves 
are from the full (full PFA) calculation. The thick solid 
curves are the full PFA limit when L/R —> cxd ; the dashed 
curves are corresponding asymptotes. 
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FIG. 15: Same as Fig. [14] but here the angle of aperture is 
kept at the fixed value 7r/2 and L/R varies. The thick solid 
curve is the full PFA result in the limit when L/R — > cxa; 
the dashed curve is the corresponding asymptote. The thin 
solid (dotted) curves are from the full (full PFA) calculation. 



In Fig. [T3]wc plot L x F {z) / tari^ {a) as function of z/L for a gold cone above a gold substrate. In PFA and full PFA 
this produces a universal curve for all sizes and for all aperture angles. We note that for large z/L the full PFA curve 
varies as {z/L)^^ which means that the force varies as z'^^ and is proportional to the projected area of the cone. For 
small relative separations the curve varies in the limit as {z/L)~^ which means that the force varies as z^^ and is 
independent of the size of the cone. In this geometry the integral equation method has convergency problems in the 
small separation limit so we omit to include any results here. This method works much better for cones with spherical 
tips which is the topic of next section. 

B. Cone with a spherical tip 

The cone with a pointed tip treated in the previous section is an idealization. In all practical cases the tip is 
rounded. In this section we let the tip be a part of a sphere. We let L be the height of the truncated cone, and R 
the radius of curvature of the spherical tip. The tip connects to the truncated cone so that the slope of the surface is 
continuous at the connection. The nearest point is z above the substrate. Let the angle of aperture be 2a. The force 
in full PFA consists of two parts, RF — RFi + RF2, where the first part comes from the spherical tip 



RFi = 27rEp (x) <^ 1 - 



and the second from the truncated cone 



X sma 



[x + 1 — sin a] 



a; + 1 — sin a 



X = z/R, 



RF2 (z) = 27rtan2 (a) Ep (x) x^ 



1 a; + 1 - sin (a) + 2L/R 

+ 1 - sin {a) (x + 1 - sin [a) + L/Rf 



z/R, 



(53) 



(54) 



where we have assumed that Eq. ((2T|) holds with 2 as an exponent. 

In Fig. [13] we compare the full PFA result (dotted curves) with that from the full calculation (thin solid curves) 
of the force on a gold cone with spherical tip above a gold substrate. In all cases L is kept 8 times the radius of 
curvature, _R, of the tip. The thick solid curves are the limiting result of the full PFA when the height of the cone goes 
toward infinity; the dashed curves are the corresponding asymptotes, RF = 27rtan^ (a) Ep (x) x. All set of curves are 
for the a values 7r/4, 7r/6 and 7r/18, respectively. In Fig. [T3]we keep the aperture angle fix with the value 7r/2 and 
vary L/R. The thin solid (dotted) curves are from the full (full PFA) calculation. The thick solid curve is the full 
PFA result in the limit when L/R — > 00; the dashed curve is the corresponding asymptote. We see that the full PFA 
results stay below the thick solid curve for all separations and comes closer the larger the L/R value. The results 
from the full non-retarded calculation cross the thick solid curve. 
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FIG. 16: Definition of the object wings. 



X. WINGS-SUBSTRATE INTERACTION 



In this section we consider a geometrical object, which we call wings, described in Fig. [111 above a substrate. The 
object has the extension L in the direction perpendicular to the plane of the figure. Thus, the bottom and top surfaces 
are squares of side length L. The thickness of the wings is 5. 



A. Thick wings 



In the standard PFA the result is independent of the thickness and is the result for infinite thickness. The potential 
is in full PFA found as 



L cos a 

V {z)^ j dSEp (w) = 2 / dxLEp {z + x tan a) 
S 

z-\-L sin a 

= 2L cot a J dwEp (w) . 



(55) 



and the force as 



F{z) = - ^^ = 2L [cot aEp {z + L sin a) - cot aEp (z)] 
= 2L cot a \Ep {z + L sin a) — Ep (z)] 



(56) 



In the non-retarded version the result is 



LF (z) = 2^2 cot aEp (z) 
= —Acqs aEp (x) ^ 
= — 2cota-Ep (x) 



{z+L sin a) 
" (l+sina/2a:) 
(1+sin a/a;)" 

(l+2a:/sina) " 

(l+a:/sin a) 



(57) 



z/L, 



where the second line is useful in finding the zero angle limit,— 4i?j, {x)lx , or the large a; asymptote, —AcosaEp {x)/x 
; the third line to find the small x asymptote, — 2cota£'p {x). In Fig. [T7]we give the results for gold wings above a 
gold substrate. The solid curves are the full PFA results for the angles 0, 5, 15, and 25 degrees. The circles are the 
corresponding results from the full integral equation calculation. In this calculation one has to choose a big, but finite 
thickness. We have chosen it to be equal to L. We have used a 2D version of the theory, i.e., we have assumed that 
the wings have an infinite extention in the direction perpendicular to the plane of the figure. We obtain the force per 
unit length and multiply this with the actual length L. One may alternatively perform the calculation for the actual 
length L directly, i.e., use a 3D version of the theory. These two results are expected to be different for large x values; 
in full PFA there is no difference. Note that the full PFA works very well for this problem. For non-zero angles the 
force is proportional to the length of the edge at small distances and to the projected surface for large distances; for 
zero angle it is proportional to the projected surface area at all distances. 
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FIG. 17: The force on thick gold wings above a gold sub- 
strate. The lines is the full PFA result and the circles is the 
full result from a 2D version of the integral equation method. 
Note that the slope of the zero angle curve is different from 
that of the others at the small x end of the figure 



FIG. 18: The force on goldwings above a gold substrate. 
The results are for angles 0° , 5° , and 25° . Thick curves (thin 
curves) are full PFA results (full results). The full results are 
for thickness 5%, 10%, 25%, 50%, and 100%. The thicker the 
curve the closer it agrees with the full PFA result for thick 
wings. The 100% curve has completely merged with the 
corresponding full PFA curve for thick wings in the whole 
zjL range covered in the figure. The full PFA results for 
thickness 5%, 10%, and 25% are indicated by circles. 



B. Thin wings 



For wings of general thickness one finds in full PFA 



F{z) 



dV(z) 



2L cot a \Ey {z + L sin a, 5) — Ep {z, 5)] 



~ 2L cot a-^ [Ep (x + sin a, 5/L) — Ep{x, 5/L)] 



z/L, 



and 



LF {x, S/L) = 2 cot a [Ep {x + sin a, S/L) - Ep {x, S/L)] . 



(58) 



(59) 



In Fig. [THlthc full PFA result for thick wings of angles 0°, 5°, and 25° is shown as thick curves. The thin curves 
represent he full result for thickness 5%, 10%, 25%, 50%, and 100%. The thicker the curve the closer it agrees with 
the full PFA result for thick wings. The 100% curve has completely merged with the corresponding full PFA curve 
for thick wings in the whole x range covered in the figure. The full PFA results for thickness 5%, 10%, and 25% are 
indicated by circles. We find that the results from the full PFA even here agree quite well with the results from the 
more well founded result from the 2D version of the integral equation method. 

If the thickness becomes extreemely small the wings behave as a bent two-dimensional metallic film, 



LF (z) ^ 2L'^ cot a [Ep {z + L sin a) ~ Ep (z)] 



« 2L^ cot a0.02A77 J nh^e'^/meVS 
= 0.04954 cot aJnh'^e'^/ me ^JJJZ 

= |a -^ 0| = 0.12385^^^^7 
= |x-.0| 
= \x^ 



{x + sin a) 



0.04954 cot a JnS.^e2/ 



oo| = 0.12385 cosaJn?i^e2/ 




{z + L sin a) ' 



-5/2 



(60) 
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FIG. 19: The force on gold wings of angles 5° and 25°, above a gold substrate at two distances as function of thickness. The 
straight lines are the full PFA result for a bent 2D metal film; the open (filled) circles are the full PFA ( 2D version of the 
integral equation method) results for the actual wing thickness, 5. 

In Fig. [in] we study how the resuhs approach the strictly two dimensional limit (straight line) when the thickness 
goes towards zero. The open circles connected with lines is the result from the full PFA. The result from the full 
calculation using the 2D version of the integral equation method is indicated by filled circles. 



XI. SUMMARY AND CONCLUSIONS 



In this work we have performed a critical test of the validity of the Proximity Force Approximation. This was 
done by comparisons with the results from calculation methods of more solid foundation, multi-pole expansions and 
a method based on an integral equation for the potentials. Calculations were performed for a variety of geometries 
including edges and corners with different orientations; we studied coated objects and free standing shells. 

We extended the traditional PFA in two ways; we did not just keep the small separation asymptote but took the 
geometrical shapes of the objects fully into account; we took the finite thickness of the coatings and free standing 
shells into account. 

Throughout this work we have presented universal figures, independent of the system scale. A word of caution is in 
place. All results are limited to the non-retarded range. This means that the range of validity in each figure depends 
on the scaling parameter. At the present status of the available measurement techniques one requires big enough 
objects, in the order of tens of micrometers, in close proximity, in the range of hundreds of nanometers. In a majority 
of cases this means that one is in the non-retarded limit and that the full PFA should be a good approximation. 

We have noted two things in particular. One is that the force between a cube and a substrate drastically changes 
character when the cube goes from standing straight to being slightly tilted. This also holds for the wings when going 
from zero angle to finite. The distance dependence of the force changes. This means that in experiments it might 
be better to keep a small but finite angle to the substrate instead of putting too much effort in trying to get perfect 
alignment; in the interpretation one should then use the full PFA result valid for the finite angle. 

The PFA results and geometrical correction factors for most of the geometries studied in this work is summarized 
in Table m 
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TABLE I: Summary of the PFA result and the geometrical correction factor for different geometries. The last column gives 
the lowest order term in the expansion of the correction factor. 



Geometry 



PFA 



Corr. Fact. 



sphere-substrate 

sphere-sphere 

oblate-substrate 

prolate-substrate 

oblate-oblate 

prolate-prolate 

cylinder-substrate 

cylinder-cylinder 

straight cube-substrate 

tilted cube-substrate 

pointed cone-substrate 

thick wings-substrate 



Rx F ^ 2nEp (x); x = z/R 

8Rx F = 2nEp (x); x = z/2R 

B-y^'^ X F = 2nEp {x) ; x = z/B 

A'y^ X F = 2nEp {x) ; x = z/A 

857"^ X F = 2nEp {x) ; x = z/2B 

SA-y'^ X F = 2nEp {x) ; x = z/2A 

R^L-^ X F= (37r/4) y/2/xEp (x) ; x = z/R 

SR^L'^ X F = (37r/4) ^j2lxEp {x); x = z/2R 

Lx F = 2Ep {x)/x; x = z/L 

L X F = Ep (a;)/(sinacosa); x = z/L 

L (tana)~^ x F = 2-KxEp {x)\ x = z/L 



Lx F 



-2 cot ctEp (x) ; 



z/L 



l/(. 


+ 1) 


l/(. 


+ 1) 


l/(a 


+ 1) 


l/(a 


+ 1) 


l/(. 


+ 1) 


l/{2 


+ 1) 

a 




a 

1 




b 


l/{x 


+ lf 


{l+2i 


/sin a) 



(l + x/sin a)'^ 



— X 

— X 

— X 
~X 

— X 

— X 

-x/A 

-x/A 



c 

-2x 
[x/sina) 



" The correction factor is: 



{'A^/tt / i)y'2TT / x{x + l)(x+2)'^ ^ x^ +2x 



'The correction factor is: 1 — cos^ a (1 + sin a/x) — sin o (1 -|- cos a/x) 
''Lowest order correction term is: — tan^ a + cot^ a\ x"^ 
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